Oscillatory transient regime in the forced dynamics of a spin torque nano-oscillator 
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We demonstrate that the transient non-autonomous dynamics of a spin torque nano-oscillator 
(STNO) under a radio- frequency (rf) driving signal is qualitatively different from the dynamics 
described by the Adler model. If the external rf current I r f is larger than a certain critical value 
I C r (determined by the STNO bias current and damping) strong oscillations of the STNO power 
and phase develop in the transient regime. The frequency of these oscillations increases with I r f as 
oc \J I r j — Icr and can reach several GHz, whereas the damping rate of the oscillations is almost 
independent of I r f. This oscillatory transient dynamics is caused by the strong STNO nonlinearity 
and should be taken into account in most STNO rf applications. 



Interest in spin torque nano-oscillators (STNOs) [l|, 0, 
H, 0] is rapidly growing among researchers as these nano- 
scale auto-oscillating systems have fascinating nonlinear 
and non-trivial properties STNOs are also interest- 
ing for rf applications due to an attractive combination 
of their properties, such as a wide range of generated 
frequencies Q, fast modulation rates 0], and easy inte- 
gration into the modern on-chip nano-electronic circuits. 

Recently, several groups have performed studies of the 
forced dynamics of an STNO under the influence of exter- 
nal microwave current [1, 0, 11 \ ■ As with other types 
of auto-oscillatory systems, such dynamics results in a 
synchronization, or injection locking, of STNO oscilla- 
tions to the external signal. The classical theory of injec- 
tion locking developed by Adler in 1946 [IH is typically 
used to analyze such experiments on STNOs. By treating 
the oscillator as an active element coupled to a resonant 
circuit, Adler obtained a simple dynamical equation for 
the phase difference tp between the auto-oscillation and 
the injected driving signal: 



^ = -Aw-Fsin(V>) . 
at 



(1) 



Here Auj = cu e — loq is the mismatch between the fre- 
quency of the injected signal u e and the frequency of 
free auto-oscillations u)q, and F is proportional to the 
amplitude of the injected signal. Despite its simplicity, 
Eq. ([I]) not only describes phase-locking phenomena in 
a variety of different physical systems, but accounts for 
all major characteristics of the phase-locking process. 
In particular, Eq. ([1]) predicts the frequency interval 
of phase-locking \Aoj\ < F and stationary phase rela- 
tion ipA = — arcsin(Aw/-F) between the driving signal 
and the locked oscillation. Additionally, it follows from 
Eq. |T]) that the phase ip approaches its locked value i/ja 
monotonically as an exponential with a time constant 
ta = 1/(F cosi/ja) inversely proportional to the driving 
signal F. 

In this Letter, we show that for sufficiently strong 
injected microwave current 7 r /, or modulation depth 



fi = Irf/Idc (where Ida is the dc bias current), Adlcr's 
model breaks down and does not give an adequate de- 
scription of the phase-locking of an STNO. The most 
striking discrepancies are i) pronounced transient oscil- 
lations of the STNO phase difference ip during its ap- 
proach to phase locking, and ii) a synchronization time 
t s which is independent of the driving amplitude I r f . We 
will show that these qualitative features are due to the 
strong nonlinearity of the STNO. Additionally, we find 
the critical modulation depth, ^ cr , separating Adlerian 
and non- Adlerian dynamics is a surprisingly small quan- 
tity. We therefore conclude that phase locking is almost 
always non- Adlerian in practical STNO devices. 

The STNO dynamics is studied numerically within 
a standard macrospin approximation [3 E ELI (H- 
The normalized (unit-length) magnetization vector m of 
the STNO free layer obeys the Landau-Lifshitz-Gilbert- 
Slonczewski (LLGS) equation QUI]: 



dm 
~~dt 



-|7|m x H e // + am x — — + |7|ajm x (m x p) . 



(2) 

Here 7 = —1.76 x 10 11 Hz/T is the gyromagnetic ra- 
tio, H e // = H a z — M s (m ■ z)z is the effective mag- 
netic field (where fioH a = 1.5 T is the external mag- 
netic field applied, in the studied case, along the normal 
z to the free STNO layer and the second term is the de- 
magnetization field with /xoMs = 0.8 T being the free 
layer saturation magnetization), a = 0.01 is the Gilbert 
damping constant, aj is the spin torque magnitude de- 
fined as aj = fn]oI /(2[ioM s eV), where K is the Planck 
constant, r/o — 0.35 is the dimensionlcss spin torque ef- 
ficiency, / is the applied current, fig is the free space 
permeability, e is the fundamental electric charge, and 
V = 3 x 10 4 nm 3 is the volume of the free layer. The unit 
vector p = cos(7o)z + sin(7o)x in Eq. ([2]), which defines 
the direction of current spin polarization, coincides with 
the magnetization direction of the fixed STNO layer. In 
our simulations we used a tilt angle 70 = 60° . 

For a constant dc current / = I^c the last term in 
Eq. ([2|) describes an effective negative damping that com- 
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pcnsatcs the natural positive magnetic damping (second 
term on the right hand side of Eq. ([2])). When the bias 
current Id c exceeds a certain threshold value I t h (2.32 
mA in our case), a stable precession of the magnetization 
vector m develops in the STNO and the free-running fre- 
quency luq of this precession depends on both H e // and 
Idc- In our simulations we used Idc = 3 mA (supercriti- 
cality parameter £ = Idc/Ith = 1-29), in which case the 
free STNO frequency was lu /(2tt) = 25.3 GHz. 

In the forced regime, when in addition to the bias cur- 
rent Idc the STNO is driven by an injected microwave 
current I r f, i.e. I(t) = Idc + Irf sin(u) e t) with ui e close to 
wo j the STNO may phase-lock to I r f. In the phase- locked 
regime the generated STNO frequency becomes exactly 
equal to ui e and a fixed (independent of the initial condi- 
tions) phase difference develops between I r t and the 
STNO oscillation. 

The results of the numerical simulations describing the 
STNO approach to phase-locking are shown in Fig. Q] for 
various I r f amplitudes. The STNO was first prepared in 
a free- running state (I r f = 0) for 50 ns to achieve a stable 
free-running regime. At t = the microwave current I r f 
was switched on with tu e = lo$. Since u> e and coq coincide, 
the phase locking manifests itself only by establishing a 
fixed phase relations between these oscillations. Fig.QJa) 
shows the time dependence of cos(0(i)) = m(t) • p (which 
is proportional to the STNO output signal) for the mod- 
ulation depth fi = I r f /Idc = 0.5. One can clearly see 
a transient beating of the envelope of the STNO sig- 
nal, which indicates an oscillatory approach to the phase 
locking. This oscillatory approach is shown explicitly in 
Fig.[l](b), where we plot the time dependence of the phase 
difference ip(t) = <fi(t) — u e t between the STNO phase 
4>{t) = arctan(TOj, / 'm x ) and the phase of the external sig- 
nal uj e t for several values of the modulation depth fx. 

We first note that the stationary value of the phase 
difference is substantially different from zero (V>o ~ 90°) 
in contrast with what one would expect from Eq. ([1]) 
for Ato = u e — ujq =0. This significant intrinsic phase 

shift [13, liz | 

is caused by the strong nonlinearity of the 
STNO generation frequency [H]. One can also sec from 
Fig.QTb) that the transient dependence ip(t) is monotonic 
only for extremely small values of fi, whereas for all the 
reasonable modulation depths strong oscillations of phase 
develop in the transient regime. The critical modulation 
depth, separating regions of the monotonic (Adlerian) 
and non-monotonic (non- Adlerian) dynamics is as small 
as fi cr = 0.0012 (curve 2 in Fig.QJb)). 

Another striking feature of the non- Adlerian transient 
regime is the very weak dependence of the synchroniza- 
tion time (the time needed to achieve the locked state) 
on the normalized driving amplitude I r f. For example, 
the envelopes of the curves 3 and 4 in Fig. QJb) have es- 
sentially the same time constant t s , whereas the classical 
Adler's model Eq. JT]) predicts that the time constant in 
these two cases should differ by a factor of 10. 
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FIG. 1: (a) Time dependence of the STNO signal cos(0(t)) = 
m(i) ■ p and (b) transient behavior of the phase difference 
between the STNO signal and I r f after I r f was switched on 
at t = 0. fi = Irf /he = 0.5 in (a) and fj, = 0.001, 0.0012, 0.01, 
and 0.1 for curves 1-4, respectively, in (b). 
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FIG. 2: Dependence on n of (a) the frequency of the transient 
STNO phase oscillations fl, and (b) its decay constant T„. 
Points - results of numerical simulations of Eq. (fSJl ; solid lines 
- analytical expressions Eq. (JS| and Eq. (0 . The vertical line 
indicates the critical modulation depth /i cr . 



The non-Adlerian STNO dynamics are further illus- 
trated in Fig. [21 where we show the dependence of the fre- 
quency of the transient phase oscillations f2 on fi (panel 
(a)) and the decay constant T s = l/r s of these oscilla- 
tions (panel (b)). One can see that Q increases with fi 
and reaches GHz values for accessible modulation depths 
fi ~ 0.5. The decay rate T s increases approximately lin- 
early with fi (following the Adler's model Eq. (JTJ) ) only 
for fi < ficr-, whereas in the non-Adlerian region fi > fi cr 
it remains virtually constant. 

The results shown in Fig. Q]and Fig. [2] clearly demon- 
strate that the STNO behavior in the transient phase- 
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locking process can not be described by a classical Adler's 
model Eq. (JTJ) . Although the transient oscillations may 
not be visible in locking experiments performed in a 
stationary regime, they may significantly influence the 
STNO operation in applications where ultra-fast transi- 
tional dynamics are important. For instance, Fig. [D^b) 
shows that the synchronization time constant t s = 1 /T s 
cannot be reduced below f=a 3 ns by increasing the ampli- 
tude of the injected current, which may limit the opera- 
tional speed of STNOs as ultra- fast signal modulators. 

To understand the origin of the non-Adlerian dynamics 
discussed above we note that STNOs, in contrast to the 
majority of conventional auto-oscillators, demonstrate a 
strong dependence of the generated frequency, lo{p), on 
the generated power p As a result, even small power 
fluctuations Sp = p — po from the free-running power 
Po may result in significant deviations of the generated 
frequency uj(j>) — u>(po) ~ NSp (here N = duj{p)/dp is 
the nonlinear frequency shift coefficient). Equations de- 
scribing STNO phase locking that also take into account 
power fluctuations have been derived in [B| and can be 
written as: 



# 

~~dt 

dSp 

~~dT 



-Alu - Fsinip + N5p, (3a) 
-2T p 5p + 2p F cos ip . (3b) 



Here ip(t) = (p(t) — uu e t is the phase difference between 
the STNO signal and the external signal, F is the nor- 
malized external signal amplitude, and T p is the damp- 
ing rate of power fluctuations. For the geometry stud- 
ied here, all the parameters of Eqs. can be calcu- 
lated analytically [f| to give: p — (£ — 1)luh / (2wm), 
N = 2u M , T p = aw H ((-l), F = ^-auj H tan(7 )/(4 % /p^), 
where £ = Idc/hh is the supercriticality parameter, 
U H = |t| (H a —M s ) is the ferromagnetic resonance (FMR) 
frequency, and u>m = \j\M s , and these expressions are 
valid for moderate supcrcriticalities £ < 1.5. For a linear 
(N = 0) oscillator Eq. (|3aj) coincides with the Adler's 
model Eq. (T]). 

The stable stationary (phase-locked) solution of 
Eqs. (J3j> has the form 



ipo — arctan(i^) — arcsin (Alo/Aloq) 



Spo = Po 



vAw + x/Awg - Aw 2 

(i + ^ 2 )r p 



(4a) 
(4b) 



where v = Npo/T p is the dimensionless nonlinearity 
parameter (in our case v f=a 1/a = 100) and Aluq = 
y/l + v 2 F is nonlinearity-enhanced frequency interval of 
phase-locking The first term in Eq. (|4aj) describes 
the above mentioned intrinsic phase shift of a strongly 
nonlinear STNO. 

By linearizing Eqs. ([3]) near the solution, Eqs. one 
can find the decay rate, A, of phase and power deviations 



from the stationary phase-locked state: 



1 



A =F P + -i^cosV'o 



±l/(r p - -Fcost/'o ) - 2^r p FsinVo 



(5) 



For a quasi- linear [y = 0), or Adlerian, auto-oscillator, 
Eq. ([5]) gives Ai = 2T p and A 2 = Fcos^o- The decay rate 
Ai describes the damping rate of the power deviations 
5p, whereas the rate A2 corresponds to the decay of the 
pure phase deviations tp — tpo . The synchronization time 
t = 1/A quantifies the overall time needed to reach a 
phase-locked state. For realistic parameters r p > F so 
the locking time of an Adlerian oscillator is determined 
by the transient phase dynamics and is given by ta = 
l/(FcosVo) oc 1/F. 

To analyze the case of a strongly nonlinear ^> 1) 
STNO, we note that in this case one can neglect F cosipo 
and simplify Eq. ([5]) to 

A = T p (l ± y/l - F/F cr ) = T p (l ± y/1 - M/Mcr) , 

(6) 

where the critical signal amplitude is F cr = 
T p /(2i/shxi()o) or, in terms of the modulation depth 
of the above considered STNO with a perpendicularly 
magnetized free layer can be expressed as 



tan7o 



(C-i) 



3/2 



2uj h 



0.0012 



(7) 



For modulation depths fi > fj, cr the decay rates A 
become complex (see Eq. (|6|)), and describe an oscilla- 
tory approach to the phase-locked state. The frequency 
of these transient oscillations is given by 



f2 = T p y/ fJ./fJ, cr - 1 



(8) 



and is shown as a solid line in Fig. H£a) for the calculated 
value T p = 0.36 ns -1 . One can see an excellent agreement 
between the numerical and analytical results. 

The decay constant T s of the phase oscillation is given 
by the smallest of the A's in Eq. (|6|) when both of them 
are real (Adlerian regime), and is equal to the real part 
of the A's, when they are complex (non-Adlerian regime): 



T p (l - fl/flcr^ , M < 



(9) 



The dependence of Eq.([9]) is shown as a solid line in 
Fig. [211b), and the agreement between the simulations 
and the analytical calculations is again remarkable. 

Thus, our analysis explains the observed non-Adlerian 
behavior of the transient STNO phase-locking as a result 
of a nonlinear coupling between the power and phase fluc- 
tuations, and provides quantitative expressions for both 
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FIG. 3: Phase-locking of an STNO to a pulsed microwave sig- 
nal with the repetition period of 16 ns. (a) Time dependence 
of the STNO signal for fi = 0.5. (b) Spectrum of the STNO 
oscillations in (a), (c) Same for /j, — 1.5. Vertical lines in (b) 
and (c) indicate the free-running frequency uq and expected 
positions of the sidebands u>o ± ^- (Other narrow sidebands 
are due to the direct frequency modulation at the pulse rep- 
etition rate.) 



the frequency Eq. ([5]) and damping rate Eq. §§§ of the 
transient phase oscillations. 

We would also like to stress that the critical modula- 
tion depth [i cri denning the boundary between the Adle- 
rian and non-Adlerian dynamics (see Eq. ((7|)), is gener- 
ally a small quantity, since a -C 1. Consequently, the 
critical microwave current I cr = fj, cr ldc is typically of the 
order of 1 /xA, which is much smaller than the typical 
injection currents used in experiments. In other words, 
phase locking of STNOs almost always takes place in the 
non-Adlerian regime. 

We finally suggest a way to observe the transient phase 
oscillations of an STNO experimentally. If the injected 
current is pulsed with the repetition rate of the order of 
l/r p , large sidebands at the frequencies ujq±Q should ap- 
pear in the spectrum of the STNO oscillations. Since the 
transient frequency f2 may be significantly larger than 
the STNO generation linewidth, both the position and 
shape of these sidebands can be measured experimen- 
tally, providing important information about such intrin- 
sic STNO parameters as the nonlinear frequency shift TV 
and the damping rate T p of power fluctuations. In Fig. [3] 
we show the results of numerical simulations of an STNO 
subjected to a pulsed rf driving signal with a repetition 
period of 16 ns. Fig. [31(a) shows the temporal profile of the 
STNO signal for a modulation depth of /i = 0.5 demon- 
strating well-resolved intrinsic oscillations of the STNO 
power. Fig. G3Jb) and (c) show the spectrum of STNO 



oscillations for two values of the modulation depth //. 
One can clearly see the sidebands caused by the intrin- 
sic transient STNO phase oscillations and their expected 
dependence on the modulation depth. 

In conclusion, we have shown that the transient forced 
dynamics of an STNO for a sufficiently strong external 
signal cannot be described by the classical Adler's model. 
The reason for this non-Adlerian behavior is the strong 
nonlincarity of the STNO generation frequency, which 
couples power and phase fluctuations. As a result, strong 
phase oscillations in the GHz frequency range develop 
during the transient regime of phase-locking. The same 
nonlinear mechanism determines the lower limit for the 
STNO synchronization time (of the order of the charac- 
teristic decay time 1/T P of the STNO power fluctuations) 
and will, likely, limit the maximum speed of the STNO 
frequency modulation. 
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